Variational principle for frozen-in vortex structures interacting with sound waves 
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General properties of conservative hydrodynamic-type models are treated from positions of the 
canonical formalism adopted for liquid continuous media, with applications to the compressible Eu- 
lerian hydrodynamics, special- and general-relativistic fluid dynamics, and two-fluid plasma model 
including the Hall-magnetohydrodynamics. A variational formulation is found for motion and in- 
teraction of frozen-in localized vortex structures and acoustic waves in a special description where 
dynamical variables are, besides the Eulerian fields of the fluid density and the potential component 
of the canonical momentum, also the shapes of frozen-in lines of the generalized vorticity. This 
variational principle can serve as a basis for approximate dynamical models with reduced number 
of degrees of freedom. 
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I. INTRODUCTION 



In many physical systems the motion of continu- 
ous liquid media can be approximately described by 
hydrodynamic-type equations having a remarkable math- 
ematical structure based on an underlying variational 
least action principle in the Lagrangian description 0, 
H H H H Q' The characteristic feature of the 
hydrodynamic-type systems is that they possess, besides 
the usual integrals of motion (the total energy, the total 
linear momentum, and the total angular momentum), 
also an infinite number of specific integrals of motion re- 
lated to the freezing-in property of canonical vorticity 
0,0,3. Thus, hydrodynamic equations describe an in- 
teraction between " soft" degrees of freedom of - frozen-in 
vortices, and "hard" degrees of freedom - acoustic modes. 
However, in the Eulerian description of flows by the den- 
sity and by the velocity (or the canonical momentum) 
fields, vortices and sound waves are "mixed". Another 
point is that due to unresolved freezing-in constraints, 
the Eulerian equations of motion do not follow directly 
from a variational principle (see for discussion). This 
work has two main purposes. The first purpose is to 
introduce a such general description of ideal flows, that 
soft and hard degrees of freedom are explicitly separated, 
and the frozen-in property of the vorticity is taken into 
account. The second purpose is to formulate a principle 
of least action in this representation. As the result, the 
acoustic waves will be described by the Eulerian fields 
of the fluid density and the potential component of the 
canonical momentum, while the canonical vorticity will 
be represented as a continuous distribution of frozen- 
in vortex lines (the so called formalism of vortex lines 
d, 0, , that previously was applied only to static den- 
sity profiles). The Lagrangian of this dynamical system is 
a non-trivial unification of the canonical Lagrangian cor- 
responding to purely potential flows, with a generalized 



Lagrangian of vortex lines. 



II. CANONICAL FORMALISM FOR FLUIDS 

A. Generalized Euler equation 

Typically in a complex classical system on the micro- 
scopic level there are permanently existing particles of 
several kinds, for instance, molecules in a gas, or the 
electrons and ions in a plasma. In general situation, dif- 
ferent components can have different macroscopically av- 
eraged velocities near a same point and/or different rel- 
ative concentrations in separated points. In such cases, 
each population of the complex fluid should be included 
into consideration individually, for example, as in the 
two-fluid plasma model discussed some later in this work. 
Now for simplicity we are going to consider the case when 
the macroscopic velocities of all components coincide and 
mutual relations between the concentrations are homoge- 
neous in space and time, so the macroscopically averaged 
physical state of the medium at a given point r = {x,y, z) 
at a given time moment t is completely determined by 
two quantities, namely by a scalar n{r,t), which is pro- 
portional to concentration of conservative particles of a 
definite sort, and by a vector j{r,t), the corresponding 
density of flow. The fleld j is related by the continuity 
equation to the field n, 



nt + divj = 0, 



(1) 



where subscript is used to denote the partial derivative. 
It is clear that j = nv, where v{r,t) is the macro- 
scopic velocity field. Let each point of the fluid medium 
be marked by a label a = (01,02,03), so the mapping 
r = x{a, t) is the full Lagrangian description of the flow. 
The less exhaustive description of the flow by the fields 
n{r,t) and j{r,t) is commonly referred as the Eulerian 
description. The relations between the Eulerian fields 
and the Lagrangian mapping are the following, 
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n{r,t) 



j 5{i — x{a, t))da, 



(2) 
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j{r,t)^ S{r - x{a,t))xt{a,t)da, 



(3) 



and they satisfy the continuity equation ^ automati- 
cally. 

With neglecting all dissipative processes (due to vis- 
cosity, diffusion, etc.), and assuming internal properties 
of the fluid homogeneous (such as the specific entropy in 
adiabatic flows, or the temperature in isothermal flows), 
the trajectories r = x{a,t) of fluid elements are deter- 
mined by the variational principle S{J Ldt) / 5x[a, t) — 0, 
with the Lagrangian of a special general form, actually 
depending only on the Eulerian fields n(r, t) and j(r, i), 

C{x{a),xt{a)} = /:{n(r),j;(r)}|„{x}j{x,xt}. (4) 

where the braces {. . .} are used to denote functional ar- 
guments as against usual scalar or vector arguments that 
are denoted by the parenthesis (...). The equation of 
motion, corresponding to the Lagrangian has a re- 
markable general structure. The usual variational Euler- 
Lagrange equation. 



d SC 



6C 



dt 5xt{a) 6x[a)^ 

in the Eulerian representation has the form (generalized 
Euler equation) 



dt 



SC 



X curl 



6C 



6C 
Sn 



(5) 



where the variational derivative SC/6j is taken at fixed 
n{r,t), while the variational derivative 5C/5n is taken 
at fixed j{r,t) (compare with Refs. [1,^, where this 
equation is written in terms of n and v and thus looks 
differently). Eq.lO) together with the continuity equation 
completely determine the time evolution of the fields 
n(r, t) and j{r, t) . 



B. Hamiltonian structure 

In the Hamiltonian description adopted for fluids as it 
is discussed in (sj, |^ , instead of the field j the variational 
derivative of the Lagrangian, 

is used (the canonical momentum). The Hamiltonian 
functional is defined as the Legendre transformation. 



n{n,p} 



6C 



j]dr~ C, 



(7) 



where j should be expressed in terms of p and n. The 
equations of motion jSJl and (0 now have the non- 
canonical Hamiltonian structure ^ , 



Pt 



nt 



- X curlp 

n\ op / 



on 



-div 



dp)' 



(8) 
(9) 



These equations can be written as p^ = {p, Ti} and 
nt — {n,7i}, where the non-canonical Poisson bracket 
is given by the following expression (see 0, 13 and refer- 
ences therein about details), 



SH 
6n 
/ curlp 
V " 





5T 


dp) 


5n 


'5T 


sn' 

Sp _ 


_ Sp 



sn 

Sp 



dr. 



dr 



(10) 



III. PARTICULAR EXAMPLES 

To emphasize universality of the employed approach, 
now we are going to consider several physically interest- 
ing examples. 

A. Eulerian hydrodynamics 

Let us start with the usual Eulerian hydrodynamics. In 
this simple case n is the density of the fluid, and the La- 
grangian is the difference between the total macroscopic 
kinetic energy and the total potential energy including 
the thermal internal energy. 



( e(n) — nUir, t) 

\ 2n 



dr. 



where s{n) is the density of the internal energy, and 
C/(r, t) is the potential of an external force. 

The canonical momentum coincides with the velocity 
field, 



p = - =v, 

n 



and the Hamiltonian is the total energy expressed in 
terms of n and p. 



n. 



p2 N 

n— + e{n) + nil {r,t)jdr. 



The equations of motion (|8I9|I with this Hamiltonian take 
the known form 



2 

Pt = [px curlp]- V(y + £'(n) + [/(r,i) 
nt = — div(np). 

B. Relativistic fluid dynamics 



In the Special Relativity (with c = 1 for simplicity) the 
field n is the time component of the (contra-variant) 4- 
vector of current, while j is the space component f^. The 
absolute value h of this 4-vector is equal to (n^ — J^Y^"^ 
and it has the meaning of the concentration of conserva- 
tive particles in the locally co-moving frame of reference. 
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The invariant expression for the action functional imphes 
the Lagrangian in the form (compare with ^^J) 



where e{n) is the relativistic density of the internal fluid 
energy including the rest energy. The canonical momen- 
tum field is defined by the relation 



(„2_ j2)l/2- 



Unlike the Eulerian hydrodynamics, now it is not pos- 
sible in general to get analytically the inverse relation 
j{n,p) in order to substitute it into the expression for 
the Hamiltonian density, 



hr = e' 



Exceptions are some special dependences e(n) (see, for 
example, @ where the particular case e{n) cx h^^^ is con- 
sidered, corresponding to the ultra-relativistic equation 
of state). 

In the General Relativity the continuity equation is 
(see Eli) 



1 



d 



ds 



where dx^/ds is the 4- velocity of the fluid element passing 
through the point (t, r), and g = det||(7ifc||, the determi- 
nant of the metric tensor gik(t,r). Therefore 



, dt 
ds ' 



dx" 
ds ' 



n = ^ {goon^ + 2ffoanj" + gc.i3rj^)/V^, (11) 

and the Lagrangian of the general-relativistic hydrody- 
namics is 



'-g.r. 



£ 
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The canonical momentum 



Pa 



e'{h) 



depends in a complicated manner on n, ji, gik- This cir- 
cumstance makes impossible in general case to present an 
analytical expression for the corresponding Hamiltonian 
functional, but, of course, it cannot cancel the existence 
of the Hamiltonian in mathematical sense. 



C. Two-fluid plasma model 

Analogously multi-component hydrodynamical models 
can be investigated where several fields n° and are 
present corresponding to different sorts of particles, with 
a = 1,2,...,^. The Hamiltonian non-canonical equa- 
tions of motion for such models have the same general 
structure as discussed above, and they should be writ- 
ten for each component. Below we consider a physically 
important example — the two-fluid plasma model. As 
special limit cases, this model contains the usual mag- 
netohydrodynamics (MHD), the Electron MHD, and the 
Hall MHD. 



1. Lagrangian formalism 

As the start point in investigation the two-fluid plasma 
model, let us consider the microscopic Lagrangian of a 
system of electrically charged classical point particles, as 
it is given in the famous book by Landau and Lifshitz |l[ . 
This Lagrangian is approximately valid up to the second 
order on v j c since excitation of the free electro- magnetic 
field by moving charges is negligible. 



2-^ Sr2 



+ 7^^^ I ^"'^^ I [va ■ Vb + {Va ■ nab){Vb ' nab)\l5) 

4^2^ \ra-rb\ V / 

a^b 

where ra{t) are the positions of the point charges Cq, 
Va{t) = ra{t) are their velocities, nab{t) are the unit 
vectors in the direction between and e;,. 



riab 



Tg - Vb 

\ra - rJ 



The first double sum in Ea. p3|l corresponds to the 
electrostatic interaction, while the second double sum 
describes the magnetic interaction via quasi-stationary 
magnetic field. It is very important that for a system 
with macroscopically huge number of particles the mag- 
netic energy can be of the same order (or even larger) 
as the macroscopic kinetic energy produced by the first 
ordinary sum in Ea. p3|l . while the terms of the fourth 
order on the velocities are often negligible. Generally 
speaking, a large part of plasma physics is governed by 
this Lagrangian, at least in the cases when the veloci- 
ties of particles are non-relativistic and the free electro- 
magnetic field is not excited significantly. Obviously, 
different physical problems need different procedures of 
macroscopic consideration of this system. The most ac- 
curate (and the most complicated) would be a kinetic 
description. However, for our purposes it is sufficient to 
apply more simple and naive procedure of the hydrody- 
namical averaging, that gives less accurate description of 
the system in terms of the concentration n{r,t) of elec- 
trons and the density j{r,t) of their flow, that satisfy 
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the continuity equation, rit + div j = (and correspond- 
ing fields N{r,t), J{r,t) for the ions, normalized to one 
elementary electric charge e, so = ZNi.) 

Neglecting all dissipative processes that take place due 
to collisions of the particles (though on this step we 
strongly reduce applicability of the following conserva- 
tive two- fluid plasma model) , we derive from Eq. (|13|) the 
following Lagrangian functional 



2/ 

" 2 



[M_j2 

.2N 
dridr2 
\ri - r2\ 
n 



+-^[curl (J-jU] , 



dr 



TpU In 



(iV(ri) - n(ri)) (^N{r2) - n(r2)) 

dr, (14) 



T-i^ In 



Z ZF(T,) 



where (J — is the divergence-free component of the 
total current. Here the constant M is the ion mass per 
one elementary electric charge, M = MijZ. The electron 
mass is TO <C M, and it will be neglected where possible. 
The magnetic energy J (B'^ /8TT)dr is included into this 
Lagrangian, where the magnetic field is 



B 



47re 



curri(j- j; 



±- 



(15) 



The terms with T^nlnn and {Ti/Z)N In N (approximate 
expressions for the densities of the thermal free energy, 
) have been introduced in order the macroscopic equa- 
tions of motion to contain the pressure terms like —'Vp/n 
(see the last term in Eq.Q), where p « nTe is the pres- 
sure of the hot electron gas, which is supposed to be 
isothermal with a temperature Tg. The functions f{Te) 
and FiTi) actually are not important since they do not 
contribute to the equations of motion and thus will be 
omitted in further equations. 

It should be kept in mind that densities of the internal 
energy should be used instead of densities of the free en- 
ergy if we suppose the flows to be isentropic. However, 
since the thermal conductivity is large at high temper- 
atures (k cx the isothermal approxi- 
mation usually works better than isentropic one. 



2. Hamiltonian formalism. Hall MHD limit 

For two-fluid models the Hamiltonian functional is de- 
fined as follows, 



V 5j 



^ 3 



SJ 



■ J j dr- £2/, 
(16) 

where j and J should be expressed in terms of the elec- 
tron and ion canonical momenta p and P. 
In our particular case we have 



P = 
P = 



MJ 47re2 2 



N c2 
mj Aire'^ 2 



curl (J - j)j 



(17) 



n 



curr^(J-j)i, (18) 



and this results in the following approximate (valid in the 
limit of small to) expression for the Hamiltonian, 



H2/ 



— {P+pf + —pj + curlp 



dridr2 



(iV(ri 



2 J J |ri -ra 

[Ten \nn+ {T,/Z)N\n N]dr 



87re2 
n[ri) I ( N{r2 



dr 

{r2)) 
(19) 



where is the potential component of the electron 
canonical momentum field. 

The equations of motion for the electron component 
now have the non-canonical Hamiltonian structure H8I9(I 
(and analogous equations for the ion component). Due to 
the large coefficient (x 1 /to in front of py in the Hamilto- 
nian, py should be small. A small value of to is also the 
reason for the electrical quasi-neutrality on sufficiently 
large scales. Therefore we can put n ^ N. Corre- 
spondingly, these simplifications will result in 7^2/ — > 
'Hsimpi{N, P + p,cmlp}. Thus, we have to deal with 
A^, with the ion velocity 

u = {P + p)/M, 

and with 

u = curlp/M = ~{e/Mc)B. 

Let us introduce the dimensionless density 

p[r,t) = N{r,t)/No, 

the ion inertial length di at the static homogeneous state, 

dl = McV(47re27Vo) = M,c^ / {'iT:Z'^e^ Nio), 

and the speed s of the ion-acoustic waves, 

s2 = (T, + ZT,)/M,. 

The Hamiltonian TY^f {p, tt, = {1/M)7isimpi for this 
simplified two-fluid plasma model (known as the Hall 
MHD) is 



Hh — 



2 1 



dr. 



(20) 



The equations of motion are 



dtu 









V 5u ) 



X curl u 



1 

+ - 
P 



, (ma 
curl — — 

SHi 



- V 



[ SHh 



dtuj — curl 



1 rSTiH ^ ^^^^^ snH_ 

p\du Su) 



X LJ 



(21) 
(22) 
(23) 
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In explicit form we easily obtain the Hall MHD equations 
as they are commonly known, 

ut = [u X curlit +— curlu; x u; -V s — -, (24) 

P 2 p 



pt = -div(pM) 



= curl 



u H — -cmluj I X u) 
P 



(25) 



(26) 



Obviously, one can normalize the length and the velocity 
scales to di and s respectively, that will result effectively 
in = I, s = 1. 

It is interesting to note that the ion generalized vortic- 
ity can be identically equal to zero (that means curl P — 
0, and therefore a; = curltt), and then we have 



Ut = 



curl curl It , 
u H X curl u 



Pt = -div(pM). 



-V^-^(27) 
2 p 

(28) 



This reduced system formally corresponds to a single- 
fluid model with the Hamiltonian 'Hs{p, u\ as follows. 



pv? (curlit) 



plnp 



dr. 



(29) 



because the equations (|27I28|) have the standard form 
similar to Eqs. H8I9(I . 



dtu 



— — — X curl u 
p\ ou J 

SHs 



Sp 



(30) 
(31) 



IV. INTERACTION BETWEEN FROZEN-IN 
VORTEX LINES AND ACOUSTIC MODES. 

The Hamiltonian non-canonical equations (|8I9|I do not 
follow directly from a variational principle. The math- 
ematical reason for this is a degeneracy of the corre- 
sponding non-canonical Poisson bracket (|10|l which is 
discussed, for instance, in 3- The degeneracy results 
in frozen-in property for the canonical vorticity field 
uj = curl p. However, representations of the canonical 
momentum in terms of auxiliary variables exist that fix 
topological structure of vortex lines, and then a varia- 
tional formulation becomes possible. A known exam- 
le of such auxiliary variables are the Clebsch variables 
ITll |. when p — V(/3 + (A/n)V/i and {n, ip), (A, p) are 
two pairs of canonically conjugate variables. But the 
Clebsch representation usually is not suitable for study- 
ing localized vortex structures like vortex filaments. Be- 
low we consider another representation for the canoni- 
cal momentum field, when dynamical variables are the 



shapes of vortex lines. For a nearly static density profile, 
n « nQ{r), such description was used in H H to study 
slow flows in spatially inhomogeneous systems. Now we 
are going to introduce a variational formulation valid for 
the general case, since the function n(r,t) is also an un- 
known variable. It will be demonstrated that variational 
principle with the Lagrangian (|41|) determines equations 
of motion for shapes of frozen-in vortex lines, for the po- 
tential component of the canonical momentum field, and 
for the density profile n{r,t). 

So, we decompose the momentum field onto the poten- 
tial component and the divergence- free component. 



p(r, t) = VLp{r, t) 



curl ^a;(r, i). 



Accordingly, the field j is decomposed, 



. 6n sn jn . 

op do; " 



(32) 



(33) 



Obviously, the continuity equation results in the relation 



(34) 



For the frozen-in vorticity field we use the so called 
vortex line representation. In the simplest form when 
the lines are closed, it reads as follows (for details and 
discussion see 0), 

a;(r, t)= d\ j> 5{r - R{u, ^ t))R^{iy, ^, t) d^ 



det\\dR/d{i^,0\\ 



R= 



(35) 



where the label v = (i^i, 1^2) G A/" belongs to a 2D man- 
ifold M and singles out an individual vortex line, while 
an arbitrary longitudinal parameter ^ determines a point 
on the line. The Jacobian of the mapping R{i', t) is 
denoted as Aet\\dR/d{v,£)\\ = {[R„i x R^^] ■ R{). 

The divergence-free component of the canonical mo- 
mentum field now is given by the expression 



, , [Re^x (r - R)]d^vd£, 
p_L = curl ^u) r,<) = / ^ ^ /, ^. 



(36) 



The vorticity variation 5u}(r, t) induced by a variation 
SR{h', S,,t) of the vortex lines takes the form Q 

6uj{r,t) = cur\r J d^v j> 5{r - R{v,S,,t))[6R x H^] d^, 

(37) 

which follows directly from Ea. (|35|l . 

It should be noted that in the case of arbitrary topol- 
ogy of the vortex lines, one has just to replace in the 
above expressions R{v,^,t) — > R{a,t), and R^d'^i/d^ 
(a;o(a) • V a)R{a,t)da, see 

Ea. H37|l results in the important relations 0], 



5n 
Jr 



(38) 
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(jjf — curlr 



Rt X fig 



dci\\dR/d{v,^)\\\ R=r 
Therefore the equation of motion for the vorticity, 

U3t = curlr [f X uj], 



(39) 



1 sn 

n op 



means 



Rt X R^ 



det\\dR/d{v, 



3\\ 



(40) 

where ^(z^i,;^2) is some arbitrary function of two vari- 
ables. A possible choice is = 0, but for general pur- 
poses we will consider below ^ 7^ 0. 

With using Eqs.||2Sl), ESI), EH), 1001), one can verify 
that if the quantities R{i',£^,t), n{r,t), and f{r,t) obey 
equations of motion corresponding to the following La- 
grangian, 

£v-s = — j niptdr — H{n, Vip + curP^u;{i?}} 

+ J S{r~R{i^,^,t)){[R^xRt]-VrA;^n)d^iyd^dr 

-J ^ivi,i'2){[Rr., X H^J • R(^)n{R)d^vdi, (41) 

then equations (|8I9|I are satisfied. Indeed, the variation 
of £v-s by SR{v, ^, t) gives the equation 



[R^ X Rt]n{R) - [R^ x j^^{R 
dR 



= JR - '^(^^^^^ 



which is easily recognized as Ea. (|40|l . 

\j\\{R)+jAR) 



Vr^(j^), (42) 



[Rt X R^] 



-hdet 



n{R) 
dR 



X Rf 



(Vr«')|r=R. 



(43) 



Variation by 5n{r, t) results in the potential compo- 
nent of the Eq.®, 



-Lpt 



Rt X Rc 



detWdR/diiy^OW 



6n 



R=r 



Finally, the variation by 6ip{r,t) gives the continuity 
equation 



nt + V • J, 



0. 



Thus, the Lagrangian (|1T|| gives a required variational 
formulation for the problem of motion and interaction 
between localized frozen-in vortex structures (described 
by the mapping R{i^,^,t)) and acoustic degrees of free- 
dom (described by the fields n(r,t) and ip{r,t)). This 



variational principle definitely can serve as a basis for fu- 
ture approximate analytical and numerical studies deal- 
ing with reduced dynamical systems where only most rel- 
evant degrees of freedom will be taken into account. 

The function ^(z^i,^2) can be useful to investigate 
nearly stationary flows, since the effective Hamiltonian, 



n = n 



Vl■V2)[[Ru^ X R^^] ■ R^)n{R)d'^v d^, 



has an extremum on stationary flows with the velocity 
field V everywhere directed along vortex surfaces. How- 
ever, one should remember that existence of globally de- 
fined vortex surfaces (and thus the function ^P) is an ex- 
ceptional case in the variety of 3D vector fields. In the 
general case one should use {u}o{a) ■ 'Va)R{o,,t)da in- 
stead of R^{i',£^,t)d^i'd^ in the Lagrangian H41|l and no 
function 5'(j/i, 1^2), since the labels ly are not defined. 

Generalization of the above theory for multi-fluid mod- 
els is straightforward. 

As an explicit example, below is given the Hamilto- 
nian of the Eulerian hydrodynamics in terms of R{v, t), 
n{r, t), and ip{r, t), 



Ti-E 



X (r - R)]d'^iydC 



dr 



(44) 



2 V ' " J 47r|r-fi|3 
+ I [e(n) + nU{r,t)]dr. 

The case of the Hall MHD 



In the Hall MHD there are two frozen-in generalized 
vorticity fields but only one density field. It is easy 
to check that in this case the corresponding Lagrangian 
takes the following form. 



>Chmhd = - y Pftdr 

+ J S{r - R{iy,^,t)){[R^ X Rt] ■VrA;^p)d^i'didr 

+ J 6{r~ Xi^i, 77, <))([X^ X Xt] ■ Vr A- V)rfV dr] dr 
-nnip, V(p -f- curr^a;{fi} -I- curr^r2{X}, u;{R}} 
-J *(i^i,i/2)([fi,.i X R^,,]- R(^)p{R)d^vd£, 

-J $(/xi,M2)([X^i X X^,]-Rr,)piX)d^pdr], (45) 

where the vector function R{i', ^, t) describes the frozen- 
in lines of the electron generalized vorticity while 
X{ii, rj, t) describes the frozen-in lines of the ion general- 
ized vorticity fi. The Hamiltonian 7i//{/3, m, a;} is given 
by Eq.|2ni). 
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